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Abstract 

We derive the free energy of the chiral Potts model by the infinite 
lattice "inversion relation" method. This method is non-rigorous in 
that it always needs appropriate analyticity assumptions. Guided by 
previous calculations based on exact finite-lattice functional relations, 
we find that in addition to the usual assumption that the free energy 
be analytic and bounded in some principal domain of the rapidity 
parameter space that includes the physical regime, we also need a 
much less obvious symmetry. We can then obtain the free energy 
by Wiener-Hopf factorization in the complex planes of appropriate 
variables. Together with the inversion relation, this symmetry relates 
the values of the free energy in all neighbouring domains to those in 
the principal domain. 

PACS: 05.50.-Fq, 02.30.Gp 

Keywords: Statistical mechanics, lattice models, chiral Potts model, 
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1 Introduction 

Onsager calculated the free energy of the two-dimensional Ising model 
by setting up an algebra that contained the row-to-row transfer matrix 
[1]. Kaufman simplified the derivation by using spinor operators (i.e. 
a Clifford algebra) [2], and Kasteleyn showed that the result could be 
obtained quite directly by pfaffians [3]. 

Most other solvable models do not appear to be amenable to such 
algebraic (in particular free-fermionic) methods. For these there are 
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three main methods of approach (cxchiding the Bethe ansatz eigenvec- 
tor method): the first and third depend expUcitly on the Yang-Baxter 
or star-triangle relation, the second implicitly. 

The first is to set up functional relations for the transfer matrix 
T. These define all the eigenvalues of T (for a finite system), so 
in principle give the exact free energy for a lattice of M rows and L 
columns. Usually one can solve these explicitly only in the large-lattice 
limit. 

The second is the inversion relation method, due to Stroganov [4]. 
In this one writes down rotation and inversion relations satisfied by 
the free energy of the infinite system. They are functional relations 
(sometimes called difference equations), wherein the free energy is 
regarded as a function of the rapidity variables (or of their difference, 
which is the spectral variable). These relations do not by themselves 
define the free energy, any more than the relation f{z + l) = 2 — 
f{z) defines a function f{z) in the complex plane. It is essential to 
make appropriate analyticity assumptions to complete the relations: 
for instance in our simple example, if we know that f{z) is analytic 
and bounded in the vertical strip < 3f?(z) < 1, then it must be 
analytic and bounded in the whole plane. By Liouville's theorem it is 
therefore a constant, so /(z) = 1. 

The third method depends on using the star-triangle relation to 
relate the free energies of the triangular and honeycomb lattices, dif- 
ferentiating this and using the fact that the derivatives are local cor- 
relations that depend only on two rapidities, rather than three. This 
gives a partial differential equation for the derivatives. As far as the 
author is aware, this method has been applied only to the Ising and 
chiral Potts models [5, 6]. 

The free energy of the chiral Potts model has been obtained only by 
the first and third methods. The result of the first [7, 8] is an explicit 
double integral, while the second [6] gives a hierarchy of equations 
that implicitly define the free energy (and can be used to obtain the 
critical exponent a) . Only recently [9] has it been verifed that the two 
results are equivalent. 

The object of this paper is to use the second method - the inversion 
relation method - to obtain the free energy of the AT-state chiral Potts 
model in the infinite-lattice limit. We do this by starting from the 
functional relations of method 1, taking the large-lattice limit, and 
showing that 

a) the resulting equations can be derived from the inversion and 
rotation relations, 

b) they can be solved by making certain analyticity assumptions. 
In fact we know these assumptions are valid from method 1, but we 
offer plausibility arguments in their favour that are independent of 
method 1. 

We then use Wiener-Hopf factorization methods to solve these 
infinite-lattice functional relations, and inevitably obtain the double 
integral result of method 1. 
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In later papers we intend to further investigate the analyticity 
properties of the free energy, or rather of the partition-function per 
site K, as a function of the two rapidity variables p and q. It is a 
mcromorphic function on a Riemann surface of infinite genus. Each 
sheet of this surface can be associated with a point on a lattice in 
2N — 1-dimensional space (for N > 2). If one fixes one of the two 
rapidities p or q, then the remaining single-rapidity space is only A^- 
dimensional. (The fact that the two-rapidity surface is only 2N — 1- 
dimensional, rather than 2N, is a manifestation of a weak "rapidity 
difference property". Unlike other planar models, it is not obvious 
that this property is particularly helpful.) 

This surface is not (for N > 2) contained in the surface generated 
by the hyperelliptic function variables known to be associated with 
the chiral Potts model [10]. Thus it is not a single valued-function 
of these particular variables, so cannot be immediately expressed in 
terms of hyperelliptic functions of them. However, the surface for the 
associated function T2{tq) is associated with a space of one dimension 
less, and is the space spanned by the hyperelliptic variables (after a 
simple transformation), so hyperelliptic functions may still be useful 
in this problem. 

There are two reasons for this investigation. One is that a large 
amount of fascinating mathematics has grown out of solvable mod- 
els, much of it connected with the explicit Jacobi elliptic function 
parametrization that is used in simpler models such as hard hexagons 
[11]. Does this generalize to the chiral Potts model and its associated 
hyperelliptic function parametrization? 

The other reason is the continuing problem of the intriguing con- 
jecture for the spontaneous magnetizations of the chiral Potts model 
[12, 13, 14]. This remains unproven, but there is an infinite-lattice 
functional relation for it [15, 16, 17] which is much like the inversion 
relation for the free energy. If one knows how to solve the latter, then 
one might hope to use similar techniques to solve the former, and 
hence verify the conjecture. Only time will tell! 

One interesting result in this paper is (71): an unexpected sim- 
ple symmetry relation satisfied by the free energy of the chiral Potts 
model. 

2 Functional relations for the transfer 
matrices 

Consider the square lattice C, drawn diagonally as in Figure 1, with 
2M rows of L sites and periodic boundary conditions. At each site 
i there is a "spin" Uj, which takes the values 0, . . . , A'^ — 1. Adjacent 
spins i and j (with i below j) interact with Boltzmann weight Wpq{ai — 
aj) for SW ^ NE edges, Wpq{ai - aj) for SE ^ NW edges. The 
parameters p, q will now be defined. 
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Let u) = e^'^*/^ be the primitive A'^th root of unity, and take k, k' 
to be two positive real constants, with A;, A;' < 1 such that k'^ + k''^ = 1. 
Also, let p = {xp,yp, Xp,tp, fip} be a set of complex parameters ("p- 
variables") related by 

kx^ = 1- k'/Xp , ky^ = l- k'Xp , (1) 

■^p ~^ Hp — ^(1 + Xp Up ) 1 tp — XpUp ) Xp = Hp . (2) 

Only one of these variables is independent. In terms of the Cp, 6p, Cp, dp 
of ref. [18], Xp = ap/dp, yp = bp/cp, fip = dp/cp. We can regard p as a 
point on an algebraic curve in (xp, y^, Ap, tp, //p)-space, and refer to it 
"rapidity" . 

Similarly, define "(j'-variables" q = {xq,yq,Xq,tq, Hq\. Then the 
Boltzmann weight functions are 

n 

Wpq{n) = Wpq{n + N) = {lJ,p/flq)''Y[{yq-UJ^Xp)/{yp-UJ^Xq) , 

i=i 

n 

Wpq(n) = Wpq{n + N) = {^ipHq)"^ JliuXp - u^Xq)/{yq - u^ypp) 

J=l 

Note that Wpq{0) = Wpq{0) = 1. We use this normalization 

throughout this paper. If Xp,Xq,yp,yq,uJXp all lie on the unit circle 
and are ordered anti-cyclically round it, then the Boltzmann weights 
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Wpq{n),Wpq{n) are real and positive for all integers n. We refer to 
this case as the physical regime. 

If Z is the usual partition function, in the physical regime we expect 
the limit 

Kj,a = lim Z^/^ (4) 

to exist, where H is the number of sites of the lattice, and the limit is 
to be taken so that the lattice is infinitely large in all directions. Wc 
refer to the partition function per site. The free energy per site 

is 

Fpq = -ksT log Kpq , (5) 

ks being Boltzmann's constant and T the temperature. Having taken 
this limit in the physical regime, wc then define Kpq outside the regime 
by analytic continuation. Wc shall find that the resulting function Kpq 
lives on an infinite Reimann surface. 

Each edge of the lattice can be regarded as associated with two 
rapidity lines, one vertical and one horizontal, as in Figure 1. The 
rapidities may be different for different lines. In 3, p is the rapidity of 
the vertical line through the edge being considered, q is the rapidity 
of the horizontal line. 

As in [19], we distinguish alternate columns of the lattice, assign- 
ing them vertical rapidities p and p' as in Figure 1. Then, as in [19], 
we can define two by row-to-row transfer matrices T, T, cor- 
responding to the two types of row of the lattice. They have elements 

L 

Ta,a' = l[Wpq{aj - aj')Wp'q{aj+i - a/) , 



Ta,a' = llWpq{aj-a/)Wp,q{aj-aj+,') . (6) 

If there are M rows of each type, then the partition function is 

Z = Trace(rf)*^ . (7) 

We regard p and p' as given, and view T, T as functions of the 
horizontal rapidity q. The parameter Hq enters the rhs of (6) only 
via its A^th power, so they are uniqely defined if both Xg and yg are 
prescribed. We write them as T{xg,yq), T{xq,yg). 

Some scalar functions of q that we shall need are 

z{tg) = yujiipiipi{tp-tg){tpi -tq)/{ypypi) ^ , (8) 

Zj{tq) = z{tq)z{utq) ■ ■ ■ z{lJ~Hq) , 
k'{l - \p\q){l - Xp'Xq) 



aq = a{Xq) 



k'>^qy^y^ 



y^y^i^p-x^) 

aq = a{l/Xq) , aqUq = ZNitq) , (9) 
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iV(Xp/ - Xq){yp> - yq){t^, - t^^ 



, _ W [ ypyp'jxp'-u^'xg) \^ 

We shall also need the by cyclic shift matrix X, with 
elements 

L 



The TT relations 

Without loss of generality, wc can take the integers k, I of [19] to be 
0,j. Then the functional relations (3.46) therein become 

T{Xq,yq)f{yq,UJ^Xg) = Tpi^qhj^q Tj (tq) + Zj{tq) TN-j{ujHq) / Uq 

(13) 

for J = 0, ... , N. 

Here Tj{tq) is the transfer matrix of the associated Tj model defined 
in eqns. (3.44) - (3.48) of [19] (which is related to the superintegrable 
chiral Potts model and to the model whose column transfer matrix is 
the Q matrix of the six- vertex model [20]). This matrix depends on q 
only via the parameter tq. in fact it is is a polynomial in tq of degree 



The T2T relation 

Similarly, equation (4.20) of [19] becomes 

{yp - U)Xq){tp' 



T2{tq)T{uXq,yq) 



UpUp' {■^p' 

LOfJ,pHp'{Xp> - LOXq)(tp 



iOtn 



T{Xq,yq) + 



XT{lj?Xq,yq) . 



I 9 \ ^ ^q, yqj . (-14) 

VpVp'iVp- ^^^Xq) 

Two other relations can immediately be obtained from (13) and 
(14) by interchanging p with p' and T with T, while leaving Tj{tq) 
unchanged. As is shown in (4.22) - (4.30) of [19], one can then deduce 
from them two sets of relations involving only the functions Tj{tq), 
which we now give. 
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The Tj relations 

The relations (4.27) of [19] axe 

Tj{tq)T2{u^-Hq) = z{J-\)XTj-l{tq) + Tj+l{tq) , (15) 
Tj{ujtq)T2{tq) = z{u)tq)XTj-i{ujHq) + Tj+iitq) , (16) 

for J = 1, . . . , A'^, where TN-^-i{tq) is defined by 

TN+l{tq) = z{tq)XTN-l{ujtq) + {aq + aq)Ti{tq) (17) 

and 

ro{tq)=0 , Tl{tq) = l , (18) 

i.e. Ti{tq) is the identity matrix. 
Invariances 

Apart from the relations k'^ + fc'^ = 1, +yp = kil + x^ Up), all the 
above equations are unchanged by multiplying Xp,Xq, yp,yq,tp,tq,k 
by e, e, e, e, e^, e^, e~^, while leaving Xp, \q, Hp, Hq and the functions 
T{xq,yq),Tj{tq), S{Xq) Unchanged. This invariancc persists if in sec- 
tion 3 we write r] and 1/r] as [(1 - k')/kf/^ and [(1 + k')/kf/^, 
respectively. 

Also, Zj{ty),aq,ag.hj q,Tj{tq),^q, S(Xq) explicitly contain the fac- 
tors /3-^, p-^,P^-\ P^-^, p-N{N-l)/2^ p-N{N-l)/2^ ^i^g^g 

/3 = {ypVpi)^ ■ These factors are "constants" (independent of q) and 
cancel out of the functional relations (13) - (17), so we could have re- 
defined the functions so as to remove them. We prefer to leave them 
in so as to make the previous invar iance more explicit. 



Comments 

Because of the star-triangle relation, the matrices T, T satisfy the 
commutation relations (2.32)- (2.33) of [19]. Prom this it follows that 
there exist invertible matrices Pi,P2, independent of the horizontal 
rapidity q, such that P{^T{xq,yq)P2, P2^f{xq,yq)Pi, P{^Tj{tq)Pi are 
all diagonal matrices, for all q. In this sense the functional relations 
above can all be simultaneously diagonalized. Their diagonal elements 
are then just scalar functional relations for each eigenvalue. 

Hereinafter we shall we shall work in this diagonal representation. 
Further, we shall only consider the eigenvalue that is the maximum 
eigenvalue of T{xq,yq)T{xq,yq) in the physical regime, so from now 
on T{xq,yq), T{xq,yq), Tjitq) are to be interpreted as the functions 
for this particular eigenvalue. From (4) and (7), noting that the the 
lattice has N = 2LM sites, the partition function per site is 



'ipq = lim 

L— »oo 



n 1/2L 

T{Xq,yq)T{Xq,yq) . (19) 
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We have distinguished the rapidities p, p' of the odd and even 
vertical rapidity hnes as it is easy to go from the above equations 
to those for a model where the rapidities p(l), . . . ,p{2L) of the 2L 
vertical rapidity lines arc all different: wherever one sees an expression 
involving p,p' raised to the power L, simply replace it by a product 
over the odd and even rapidities, respectively. Thus (11) becomes 



However, from now on we shall only consider the fully homogeneous 
model, with vertical rapidity p' = p for all columns. In this case T 
and T differ only in a cyclic shift of the L spins, and since we are 
considering only the maximum eigenvalue, the associated eigenvector 
is unchanged by such a shift. This eigenvector is also unchanged by 
multiplication by X, so hereinafter we take 



3 Functional relations in the infinite- 
lattice limit 



We emphasize that the above equations are exact for a lattice with 
a finite number L of sites per row. In [7, 8, 21] they were solved in 
the limit of L large. A key step in this working was to note that in a 
selected domain on the {xq,yq) surface, certain terms in each relation 
became exponentially small (relative to the other terms) as L became 
large, so could be neglected. 

We can identify such terms by examining the low-temperature 
limit, when k' ^ and k ^ 1. Similarly to [7], take yp,yq — >■ 1, 



Xp ~ tp, Xq ~ tq, Xp = 0{k'), Xq = 0(fc')- Then z{tq) = 0{k'^^'^), 



T{Xq,yq) ~ 1 , T2itq) ~ (1 - C^t,)^ , Tjitq) ~ 0(1) . (21) 



It follows that all terms in (15) - (17) that contain the function 
z{tq) or ctq will be negligible, giving 




p'=P , f{Xq,yq)=T{Xq,yq) , X = l . 



(20) 



and from the definitions in [19] 



Tj{tq) = T2{tq)T{utq) ■ ■ ■ T2{U^ \) 



(22) 



T2{tq)T2{ujtq) ■ ■ ■ T2{u;'^ ^tg) = 

Similarly, the last term in (14) will be neglible, so 



(23) 



r2{tq) 




T{Xq,yq) 
T{uXq,yq) 



(24) 
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Replacing Xg,yg,tq in this equation by u;^ Xg,yq,LoHq and taking the 

product over j = 0, . . . , A'" — 1, we obtain (23). 

The last term on the rhs of (13) will also be neglible, provided 
= o{k'^^~^^''^^), i.e. for j = 1, . . . ,N. (The j = equation is the 

same as j = N, so wc lose nothing by ignoring j = 0.) 
Taking j = 1, we obtain 

T{Xq, Vq) T{yq,UJXq) = Vp^q = rp{Xq, yq) . (25) 

Keeping yq fixed, replacing j by ,7 — 1 and Xq by ujXq in (13), then 
dividing the resulting equation into the original equation (13), we get 

, ^^^^^"/'^^ X r^itq) (26) 

{yp-U:Xq){tp-tq)\ ^ ' 

for j = 2, . . . ,N. But this is the same as (24), so there is only one 
further equation contained in the set (13), namely (25).^ 

To summarize: (24) can be regarded as defining T2{tq), and (22) as 
defining Tj{tq). The product relation (23) is a direct consequence of 
these definitions. So, out of all the functional relations we originally 
wrote down, the only relation left that contains information about the 
function T{xq,yq) is (25). 



T{Xq,yq) 

T{ivxq,yg) 



The function S{\g) 

We shall need another function, also introduced in [7],^ namely 



SiXq) = ^^T{Xq,yq)T{LJXq,yq)---T{u''-^Xq,yq) , (27) 



where 



JV-l 

= n 



(28) 



In [7, 8] it is shown, for finite L, that ^(A^) is a polynomial in Xq of 
degree {N — 1)L. 



Domains 

The parameters Xq,yq,tq, Xq are multi- valued functions of one another 
and we have to be careful to identify them when working with the 
above equations. 

Prom (1) and (2), tq and A^ are related by 

kh^ = l-k'{Xq + Xq^) + k''^ . (29) 

^To put this point another way, in the large lattice limit all the N + 1 identities (13) 
can be obtained from (24) and (25). 

^This definition differs from that in [7], but only in an extra factor N^/'^ 

(_l)(Af-2)(JV-l)/2^ 
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Figure 2: The complex tg-plane with its N branch cuts Mq, . . . , M/v-i 

(for iV = 3). 

Define 

r? = [(l-A;')/(l + A;')]'/'^ , < r/ < 1 , (30) 

and consider the cut t^-plane shown in Figure 2, with branch points 
at tq = uj^r] and tq = lo^ /i] and a brancli cut Mj in between, for 
j = 0, . . . N — l,as shown. Then Xq is a single- valued analytic function 
of tq in this cut plane, and we can require (consistently with the low- 
temperature regime considered above) that |Ag| < 1. We can then 
choose yq to be the solution of (1) such that | argy^l < tt/2N. (If k' 
is small, this means that yq ~ 1.) Finally, we define Xq to be tq/yq. 

The result is that Xq lies in the region £ of Figure 3, while yq lies in 
the region TZq. The boundaries of TZq, . . . ,TZn-i arc where \\q\ = 1. 

Let us use the symbol Pq to denote the domain we have just spec- 
ified, i.e. 

Vq : \\q\ < I , Xq e S , yq e TZq . (31) 

Since tq and Xq are uniquely determined if both Xq and yq are known, 
we shall often say that ^\xq, yq) lies in Vo" if the above constraints are 
satisfied, but it must of course be remembered that Xq , yq are related 
complex variables. 

If k' is small, /c ~ 1 and {xq,yq) lies in Vq, then yg 1 and Xq is 
barred only from small regions about the N points l,uj,. . . u!^~^. 

If Xq is an allowed value in Vq, then so are co^Xq, for all integers j. 
We shall sometimes regard Xq as the independent variable in Dq. 

We shall need to consider the domains (Riemann sheets in the tq- 
variable) neighbouring Vq. There are N of these, obtained from Vq 
by moving tq across one of the N branch cuts in Figure 2. Then Xq 
moves into one of TZq, . . . ,TZn-i in Figure 3, while yq moves into 8. 
For instance, if we move tq across the cut from uj'ri to w'^'/ry, we 



10 




Figure 3: The + 1 regions £,TZo, . . . ,TZn-i of the complex plane 
in which Xg and t/g lie (for = 3). TZq, . . . TZn-i are the 
interiors of the approximate circles centred on 1, a;, ... , 
E is the the complex plane outside all N circles. 

enter the domain P^, where 

: lA^I >\ , Xg^Tlr , Vg^S . (32) 

Now it is Xg that is severely restricted,, while is free to move about 
most of the complex plane. 

By continuation from the low-temperature, small A;', limit, the 
functional relations (22) - (28) hold for {xq,yg) and (xp,yp) both lying 
in Vq, except that we must be particularly careful with (25), since 
the arguments {xg,yg), {Xg,y'g) = {yg,ujXg) of the two T-functions 
therein cannot both lie in Dq. One of them must lie in a neighbouring 
domain, but which domain? The answer is that we want both (xg, yg) 
and {x'g,y'g) to be in or near the physical regime. Only there can we 
expect to obtain the correct free energy in the large-lattice limit. 

If {xq,yq) is in the physical regime, then, from our remarks before 
(4), Xq,yg,tq lie on the unit circle and 

argyg - 27r/A/" < arga;^ < argy^ , 

If it also in Vq, then yq must be inside TZq, so Xq must be on that 
portion of the unit circle between between regions TZq and TZn-i in 
Figure 3. Hence Xg G TZq, so {x'q,y'q) lies in Vq. Writing Tr{xq,yq) 
for the analytic continuation of T{xg,yq) to domain V'^, we should 
therefore write (25) as 

T{Xq,yq)To{yq,UJXq) = rp{xq,yq) . (33) 

Alternatively, {x'q,y'q) may lie in the physical regime and in T)q, 
with ?/g G TZq. Then Xq G TZn-i and, dropping the sufHxes on (x^, y^), 
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(25) becomes 

TN-l{uJ~'^yq,Xg)T{Xq,yq) = rp{uj~'^yq,Xq) . (34) 

Both (33) and (34), with {xq,yq) G Vq, are valid interpretations of 
(25). 

We refer to T>q as the principal or central domain (or Riemann 
sheet). In a subsequent paper we intend to discuss the full Riemann 
surface formed by analytically continuing lim^i^^oo T{xq, yqY^^ in both 
the q and the p variables. 



The inversion relation 

We now turn from the first to the second general method for calculat- 
ing free energies, the inversion relation method [4, 22, 23, 24]. Let 

N-l 

Dpg = det^Wpqii - j) , Ppq=l[ Wpq{n) , (35) 

n=0 

i.e. Dpq is the determinant of the N hy N Toeplitz (cyclic) matrix 
with entry Wpq{i—j) in row i and column j. Then under quite general 
circumstances the partition function per site Kpq satisfies the inversion 
relation 

-il/N 



D P D P 

PQ pq-'-^qp-'- qp 



(36) 



For the chiral Potts model PpqPqp = 1 and (36) becomes, using 
(2.48) of [19], 

^pql^qp — fpq i (^''^) 

where rpq is defined in (10). 

Let R be the operator that acts on the rapidity p so that 

XRp = yp , yRp = ujxp , HRp = l/np . (38) 

Then replacing j*, g by Rp is equivalent to rotating the lattice through 
90°. This does not change the partition function per site, so 

l^pq — l^q,Rp • (^9) 

Combining (37) and (39) gives 

Kpqi^p,Rq = ■ (40) 
Remembering that Kpq = T{xq,yq)'^/^, it follows that 

T{Xq,yq)T{yg,UJXq) =rpq . (41) 

This is precisely the relation (25), again with one function T be- 
ing understood to be the analytic continuation of the other from the 
physical regime. Thus we did not need to go through all the working 
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that was necessary in [20], [19] to derive the functional relations (13) - 
(18). We could have used the simple inversion relation, obtained (41), 
i.e. (25), then defined the auxiliary functions Tj{tq), S{Xq) so as to 
obtain the infinite-lattice relations (22) - (28) above. 

4 Solution of the infinite-lattice func- 
tional relations 

We now solve the infinite-lattice relations (22) - (28) successively for 
'^2(^9)) S{Xq), T{xq,yq), thereby obtaining the partition function per 
site Kpq. Each side of every one of these relations is a quantity raised 
to the power L, and this is the only way L enters the relations. We 
should define new functions that are the Lth roots of T2(ig), S{Xq), 
T{xq, Uq), and write down the Lth roots of each relation. Equivalently, 
from now on we formally set 

L = 1 (42) 

in equations (8) - (11) and (21) - (41). 

Thus T{xq,yq), or Tpq, is now the partition function per site Kpq 
of the lattice, defined as in the paragraph containing (4), (5), i.e. 

Kpq = Tpq = T{Xq, yq) . (43) 

Our derivation will be based solely on the inversion relation (25), 
as interpreted in (33), and the definitions in (22) - (28). It will nec- 
essarily involve various analyticity assumptions, and in making these 
we have been heavily guided by the previous derivations in [7] , [8] . In 
particular, we shall regard T2{tq) as a function of tq (rather than, say, 
Xq or ^g), and S{\q) as a function of Ay. It must be admitted that 
from the present point of view this is by no means an obvious thing to 
do. We shall present what justifications we can, and attempt to state 
clearly the analyticity assumptions that we make. 

In particular, the domains Dqi • • • ' ^7V-i ^'"^ made up of + 1 
connected Riemann sheets of the cut tq plane of Figure 2. Together 
they are the beginnings of the full Riemann surface on which the 
functions live (and which we intend to discuss in a later paper). Here 
we only need these domains or sheets. In fact we only need the surface 
consisting of Vq and the adjacent neighbourhoods of V'q, . . . ,'D'j^_^ 
(obtained by just crossing the branch cuts in Figure 2). Since |Ap| < 1, 
and on Vq it is true that \Xq\ < 1, we can choose this surface so that 

[ApAgl < 1 on P+ . (44) 

Let us call this extended surface P+. It includes the central domain 
Pq and all functions are defined on it by analytic continuation from 
Vq. Then a basic assumption that we make is: 
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ASSUMPTION 1: The function T{xq,yq) is non-zero and analytic 

on 



This is consistent with the low-temperature result and is a stan- 
dard assumption in the inversion relation method: the free energy 
is analytic in some fundamental domain that includes the physical 
regime. The remark after equation (7) that fig enters T{xq,yq) only 
via its iVth power Xg is important here: it it were not true then there 
would also be branch cuts at and oo in the complex tg plane, and 
the surface would be even more complicated. 

In [7, 8] wc took L to be large but finite and used the fact that 
Tj{tq), S{Xq) are then polynomials to solve the functional relations. 
Here we adopt a different (but related) strategy: we use Weiner-Hopf 
factorization to solve (22) - (28). 



Calculation of T2{tg) 

Since P+ extends beyond the branch cuts in Fig. 2, T2{tq) is not 
necessarily a single- valued function of tq-. we may expect it to have 
these branch cuts in the complex tq -plane. 

First consider possible poles or zeros. From Assumption 1 and (24), 
these can occur only when Xq = \p, Xq = Xp and yg = y^, or when 
Xq = 1/Xp, Xg = yp and yg = Xp . The second possibility is excluded 
by (44), so we have to consider the first. Since yp,yg G TZq, this can 
only happen when yg = yg. But then the term {tp — tg)/{xp — Xg) 
simplifies to yp: the potential pole and zero cancel one another. Thus 
T2(tg) is analytic and non-zero in From (22), so therefore is rAr(tg). 

Prom (22) and (24), 



TN{tq 



{xp-uj ^Xg){y^ ~ x^){t^ -t^) T{xg,yg 



' yp^'-^ix^ -x^){yp-Xq){tp-uj-Hq)T{uj-^Xg,yq) 

(45) 

Also, from (10) and (33), 

1 ^ « - x^){y^ - y^){tp - co-Hq) Uyq,xq) 

T{u-^xq,yg) N{xp-u-^xq){yp-yq){t^ -t^) 

(46) 

Combining these together, there are several cancellations, leaving 



^^^^ ^ ^ {Vp - Xq){yp -y^)T{xg,yg) Tojyq, x^, ^^^^ 
^ N yf^'"^ {yp- Xg){yp- yq" 



Now To(xq,yq) is the analytic continuation of T{xq,yg) across the 
cut Mo, where Xq and yg are both on the boundary of TZq. Analyti- 
cally continuing the rhs of (47) across this cut and then interchanging 
Xq with yg leaves it unchanged, while also leaving tg unchanged. For 
a given tq, these are the only possible values of {xq,yq) in this neigh- 
bourhood. Thus in the neighbourhood of Mq, on either side of the cut, 
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TNitq) is a single- valued function of tg. There is therefore no need for 
the cut Mo for this function: it can be removed. 
Prom (22), the equation (23) can be written as 

TN{tg)T2{uJ-'^tg) = aq . (48) 

The rhs of this equation is certainly is not a single- valued function of 
tq across the cut Mq (when Xq l/Ag), so nor is T2{io^^tq). Hence 
the plane in which T2{tq) is analytic must have M;v-i as a branch cut. 

ASSUMPTION 2: T2{tq) is single-valued, with logarithmic deriva- 
tive zero at infinity, in the tq plane containing only the cut Mjv-i . 

This is by no means obvious. It is consistent with Ti\f{tq) = 
T2{tq) • • • T2{oJ^^'^tq) uot liaviug the cut Mq, but is not implied by 
it. It is possible that T2{tq) could have cuts other than Mjv-i that 
cancelled one another in the product function TN{tq). On the other 
hand, it is the simplest assumption consistent with T]\f{tq) not having 
Mq as a branch cut, and again part of the spirit of the inversion 
relation method is to assume the simplest possible analytic struc- 
ture.^ Exhibiting the possible multi-valuedness of T2{tq) by writing 
it as T2{xq,yq), our assumption implies that for {xq,yq) on I>+, 

T2{Xq,yq) = T2{uj'' yq,Uj"'' Xq) (49) 

for r = 0, . . . , — 2. So we are assuming these additional symmetries. 

The requirement that the logarithmic derivative be zero at infinity 
follows from the finitcness and analyticity of the Boltzmann weights 
and hence T{xq, yq) when Xq ^ co va. £. 

It follows from these two assumptions and the above remarks that 
there is a closed curve C is the complex tq plane, surrounding the 
potential branch points r/, l/ry as in Figure 2, such that 

i) T]sf{tq) is analytic and non-zero inside and on C, 

ii) T2{uj~^tq) is analytic and non-zero outside and on C, 

iii) {d/dtq) \ogT2{oj~^tq) — *■ as — ^ GO. 

We can now solve (48) by Wiener-Hopf factorization. [25] Regard 
Uq as a function a.{tq) of tq and temporarily drop the suffix q. Define 

F_(s) = --^ / ^— f-f-loga(t)) , s outside C , (50) 
2m Jc t — s \dt J 

F^{s) = ^/^^(|log«w)cit , s inside C . 

Shifting the contour for -F-(s) to be inside C, that for -F+(s) to 
be outside, and s inbetween, it follows from Cauchy's integral formula 
that 

F_(s)+F+(.s) = a'{s)/a{s) . (51) 

^We of course know that Assumption 2 is true from [7], but are trying to present an 
argument based on the inversion relation. 
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Hence from (48), 



F_{s) - ^logT2{u;-h) = -F+is) + -^logriv(s) . (52) 
as as 

The rhs of this equation is analytic for s inside and on C, while the 
Ihs is analytic outside and on C and tends to zero as s ^ cx). Hence 
both sides are entire and bounded. By Liouville's theorem they must 
each be constant, in fact zero, so 



^l0gT2{u-'s)=F^{s) 

as 



(53) 



We can shrink C to just surround the branch cut in the t plane 
from 77 to 1/r]. Then as t goes round C, A = Ag (as defined by (29), 
goes anti-clockwise round the unit circle. Changing the variable of 
integration in (50) to 6, where A = e*^, then replacing s by ujtq, it 
follows that 



logr2(t,) = ^ 



•2'^ / 1 + XpcJ 



1 - Xpe'^ 



log 



A{e) - iOtq 



de 



(54) 



for tq lying in the complex plane with a single branch cut from uj ^rj 
to /rj. Here 



'1 -2k' cose + k'^Y^^ 



fc2 



(55) 



and we have integrated (53). We have fixed the integration constant 
by using (23) and the formula 



2tt Jo \l 



J8" 



iog[Aief-t''] = log 



k'jl - XpXqY^ 

k'^Xn 



(56) 

which is true when tq, Xq are related by (29) and |Ap| < 1, \Xq\ < 1. It 
also follows from this formula and (48) that TN{tq) = T2{tq) ■ ■ ■ T2{u}^~'^tq), 
in agreement with (22). 



Calculation of ^'(A^) 

The function S{Xq) is defined by (27) in the domain Dq, where |Aq| 
is sufficiently small as to justify our neglect of various terms in the 
functional relations in the limit of L large. (At low temperatures, 
this means that |Aq| < k'). Outside this domain we here define it by 
analytic continuation. 

This is a different definition from [7], where much use was made 
of the fact that S{Xq) is a polynomial for finite L, with zeros located 
approximately on concentric circles between |Ag| = k' and \\q\ = l/k' . 
In the large-L limit these cause the function S{Xq) of [7] to have a 
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different analytic form across each such circle of zeros, so the two 
definitions will only agree for \Xq\ inside the smallest circle. 

We first repeatedly use (24) to express T{ujXq, yq), . . . , T{uj'^~^Xq, Uq) 
in (27) in terms of T{xq,yq), giving 

(57) 

Within Vq, if Xq = fi^ is given, we take Uq to be the root of (1) lying 
in TZq. There arc N possible choices for Xq in £, but by the construction 
of (27), S{Xq) is the same for each, so it is a single- valued function of 
Xq in Vq and (by analytic continuation) in P+. Prom Assumption 1 
and the definition (27), it is analytic and non-zero in P+. 

For |Aq| close to one, choose Xq, like yq, to lie near the boundary 
of TZq. Replacing by l/Xq simply interchanges Xq with yq, while 
leaving tq unchanged. Note that T2{uj''^tq) does not occur on the rhs 
of (57), so allowing Xq to move just outside the unit circle docs not 
take any of the T2 functions in (57) across a branch cut in the tq plane 
and they are unchanged by analytically continuing from to 1/Xq. 

The function T{xq,yq) is simply replaced by To{yq,Xq), which is 
given by (47). So analytically continuing (57) from A^ to 1/Xq (thereby 
interchanging Xq with yq), then multiplying by the original equation, 
we obtain (after many cancellations) 

S{Xq)S{l/Xq) = (A,/,f ) n \^^j-hJ)N-2j ■ (58) 

This is not the same as equation (26) of [7], because here S{l/Xq) is 
defined by analytic continuation of S{l/Xq) from |A| < 1 to |A| > 1. 
As remarked above, this is a different definition from that used in [7]. 

The rhs is known. As Xq moves round the unit circle, tq moves on 
the positive real axis from i] to l/rj and back again. The logarithm of 
the rhs is analytic on this contour, returning to its original value. It 
follows from the above remarks that log S{Xq) is analytic inside and 
on the unit circle in the complex Xq-plane. 

We can therefore solve (58) for S{Xq) by Wiener-Hopf factorization. 
Writing the rhs of (58) as R[tq], we obtain 

where A{9) is defined by (55). Define,^ for |Ap| < 1 and \Xq\ < 1: 



*These definitions are the same as those of eqns. (45) and (46) of [7], except for the 
factor k"^^^ — inside the logarithm in Apq. This gives an extra constant contribution 
(AT - 1) log /c-h iTT{N - 2){N - 1)/12 - {N/2) log AT to Apq, and ensures that Apg = in 
the low-temperature limit of the Appendix. 
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B. 



pq 



1 

8^2 Jo 

N-l 



2lT f2TT 



1 + A„e^^ 1 + A„e^ 



1 - Ape*^ 1 - AqC* 



X 



^ (AT - 2j) log[a;-^/2A(^) - J/'^A{(t))] dO d(j) , (61) 



it follows that 



log5(Ag) = (iV-l)log 



■^qp Bpq 



(62) 



Let q = Xq, 1/Ag, ig, l/zig}, so Ag = 1/Ag and tq = tq. One can 
verify that (62) does indeed satisfy (58) by analytically continuing 
Aqp,Bpq to |Ag| > 1 (deforming the contours of integration accord- 
ingly), and verifying that the resulting functions satisfy 



-^qp "I" -^'i 



-,P - 2^(^^-j) log 

N-l 



(63) 
(64) 



Bpq + Bp-q = 5^ (AT - 2j) log 

i=i 

provided that for A^ near (on) the unit circle, tq is taken to be the 
solution of (29) near (on) the real positive axis, so that |arg(iq)| < 
tt/N. 



Calculation of T{xq,yq) 

Now wc know S{Xq), we can solve (57) for Hpq = T{xg,yq). The result 
is best expressed in terms of T{xq^yq)/{DpqPpqY/^ , where Dpq, Ppq 
are defined in (35). Prom (2.44) of [19], 



N-l 



Dpq = n 



{1 - iO^){tp - CvHq) 

{xp -u;ixq){yp-u^yq) 



while from (3), 



P. - (Ap/A,)(^-^)/^ n 



yq - LQJXp 

yp - uixq 



Using these formulae and (28), (57) gives 

TiXq.yq)'' 



D P 

pq pq 



^-1 f 1-1+ \ 

(Af«,)-^5(A,)n 



N-j 



= S{Xq) exp(A 



{Xq/Xp 



PIJ 



,(7V-l)/2 



k'{l — ApAq)2 



N-l 



eyip{ Apq — Aqp — Bpq) 



(65) 



(66) 



(67) 
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This is the result (50) of [7]. 

There are two "inversion" relations: the inversion relation (37) and 
the rotation symmetry (39). We have derived (67) from the single 
combined relation (41). However, Bpq = —Bgp, so the rhs of (67) is 
inverted by interchanging p with q. Thus (37) is indeed satisfied and 
we have in fact satisfied both relations. 



5 Analytic continuation of T(xg, i/q) 

We can readily use the above equations to obtain the analytic continu- 
ation of T{xq, Uq) from Vq to the neighbouring domains V'q, . . . , V'j^_^. 
As in section 3, if {xq, yq) G Pq we write T{xq^ yq) simply as T{xq, yg), 
while if {xq,yq) G i^*^^ — 0, . . . ,N — 1) we write T{xq,yq) as 

Tj.[Xq, yq)- 

First note by repeated use of (24) that, for {xq,yq) in Vq and 
r = 0,...,N-l, 

T{u^Xq,yq) = T{Xq,yq) f[ " ^'"l^^ " ^ " ^^^^ 

Now analytically continue {xq,yq) from Pq to Pq) that tq goes 
through the branch cut Mq. Then {uj'^Xq, yq) on the rhs goes from Pq 
through the cut Mr to P^. The function T{xq, yq) on the rhs therefore 
becomes the function To{xq, yq), while T{u'^Xq, yq) on the rhs becomes 

Tr{uj''Xq,yq). 

Going through the branch cut Mq and returning to the original 
value of tq is equivalent to interchanging Xq with yq. (The branch 
points of Mq are where Xq = yq.) We can therefore interchange Xq 
with yq and obtain the equation 

Like (68). this equation holds for {xq,yq) in Pq and r = 0, . . . , N — 
1. The functions T2{tq), . . . ,T2{uj^'~^tq) in these equations are single- 
valued and analytic functions of tq across the cut Mq, i.e. Mq is not a 
brach cut of these functions. Hence they are the same functions, with 
the same values, in each equation. 

We can therefore eliminate the T2 functions by dividing (69) by 
(68). Using also the form (33) of the inversion relation , we obtain, 
for (xq, yq) in Pq and r = 0, . . . , A?^ — 1, 

TUV X) = ^Tiu^^q,yq) - tp-U^-Hq 

n yg, g) T{oj-^Xq,yq)T{xq,yq) jj-^{xp-u;^-%){yp-ujxq) 

N-l . 

X n - — . ^ . (70) 

[Xp-OJ^ ^Xq)[yp-UJ^yq) 



i=r+l 
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Given T{xg, yq) in Dq, this relation enables us to calulate the func- 
tion in each neighbouring domain P^. If we set r = we regain 
the inversion relation (33), while if we set r = A?" — 1 we obtain the 
alternative form (34) predicted in section 3. 

One can continue. If we allow Xq in (70) to move from £ into TZg 
(with s = 1, . . . , A'' — 1), then {x'q,y'q) = {uj'^yq,Xq) on the rhs moves 
into a new domain Vrs, where x'qES and y'q £ TZg- Correspondingly, 
Tr{x'q,y'q) On the Ihs becomes the function Trs{x'q,y'q). On the rhs, 
the T functions become Tr+s,Ts-i,Ts, all of which can in turn be 
expressed in terms of the original function T (in domain Pq) by using 
(70). 

The wc can continue x'q into TZt to obtain the function Trst whose 
arguments lie in the domain T>j.g^, and so on. If each function were 
unique, then this procedure would give a Cayley tree of functions and 
domains, each having N neighbours, with no circuits. This would 
form an infinite-dimensional Bethe lattice. In a subsequent paper we 
intend to show that in fact the functions are not unique and there 
are circuits, so that one obtains instead an A^'-dimensional lattice (for 
N > 2), each site corresponding to a domain and a function. This is 
the Riemann surface on which T{xq, yq) lives. 

By contrast, the Boltzmann weights and related quantities of course 
form a zero-dimensional surface: for instance, the function Dpq in (65) 
has just 2N sheets in the tq plane, corresponding to either Xq or yq 
being in TZr (for r = 0, . . . , A'^ — 1), the other in £. 

We also hope to explore the question of whether one can describe 
the Riemann surface for T{xq, yq) by an appropriate generalization of 
the hyperelliptic variables and functions defined in [10]. 



6 Analytic continuation in Xp, yp 

Up till now we have held the p variables (corresponding to the vertical 
rapidity) fixed at some values in Vq, except possibly in the inversion 
relations (36) - (38). One can of course vary the p variables as well as 
the q variables, and this would be necessary for a full understanding 
of the Riemann surface on which then free energy lives. In particular, 
we should like to obtain the variable analogue of (70). 

Prom this point of view T is a function T{xp, yp\xq, yq) of both the 
p and q variables. It has an unexpected ratio property. From (54) it is 
apparent that T2{tq) is unchanged by replacing Xp,yp by u;Xp,yp (the 
variables remain in Vq and Xp is unaltered). Prom (24) it follows that, 
for ixp,yp) and {xq,yq) both in Vq, 

TjXp, yp\Xq, yq) T{uJXp, yp\u>Xq, yq) ^ {Xp - Xq){ujtp - tq) 
T{Xp,yp\uJXq,yq)T{u;Xp,yp\Xq,yq) {uJXp - Xq){tp - tq) 

This result is not at all obvious from first principles: it could serve as 
a useful test of series expansions. 
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Note that one has to be very careful in applying the equations 
of this paper in domains other than those for which they were ob- 
tained. This relation is a good example: if we analytically continue 
T{xp,yp\xg,yg) to {xp,yp) £ Vq, {xg, yg) G P'. by allowing Xg to move 
from £ to TZr, yq from TZq to then the T- functions with third ar- 
gument Xq will become the functions above, while those with third 
argument uxq will become T^+i. The different functions, originally on 
the same sheet, will move onto different sheets. 

For the rest of this section wc shall hold a;^, yq fixed, lying in Pq; 
so for brevity we write T{xp, yp] Xq, yq) as T[xp, yp]. 

Fortunately it is not necessary to repeat all the above working. 
We can replace Xp, yp in the final result (67) by wxp, yp (so again the 
variables Xp, yp remain in Dq) and take the ratio of the equation with 
these arguments to that with the original arguments. Since Xp is the 
same in both, the factors involving Apq,Bpq cancel, and we obtain 

^ /. X jtq - u;tp){yg - uJXp) T[u)Xp,yp] 

^2(tpj = 2? ^ ^[ r ' y'^) 

VqiXq-UJXp) T\:zp,yp\ 

where f2(tp) is the function T2{tq) defined by (54), but with p inter- 
changed with q, i.e. 




A(6l) - uotp 

Vq 



de . (73) 



It is analytic in the tp plane, except for a single branch cut from from 
oj~^ri to oj~^rj~^. This result is of course consistent with (71). 

Let Tr[xp, yp] be the analytic continuation of r[a;p, yp] from (xp, yp) G 
Pq to (xp, yp) G V'^ (i.e. Xp moves to the region TZ^- of Figure 3, yp to 
the region £). Then from the inversion relations (36) - (38), 

T[u)''^yp,Xp]TQ[xp,yp] = Tpq = TN-i[u)''^yp,Xp]T[xp,tjp] . (74) 

(The second of these relations is the analytic contination of the first.) 

Proceeding as in the previous section, from (72) and either equa- 
tion of (74), we can deduce that 



NT[u''xp,yp] tg - LoH^ 



p 



T[uj ^Xp,yp]T[xp,yp] jj-^ {xq-ixi^Xp){yq-uj^yp) 



N-l , j. 

X n - — . ^ (75) 

j=r+l K ~ ^^VpAVq - ^^Xp) 

for r = 0, ...,Ar- 1. 

If r = we regain the first inversion relation (74); if r = AT — 1 we 
obtain the second. 
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7 Summary 



A basic difficulty with the A''-state chiral Potts model is that it lacks 
the "rapidity-difference" property [8]. Because of this, there is no 
known explicit parametrization in terms of elliptic or other single- 
valued single- variable functions. The author did show that one could 
parametrize the Boltzmann weights in terms of hyperelliptic functions 
with N — 1 related arguments [10], but so far little progress has been 
made using this parametrization. (For A'^ = 3 it has been used to 
obtain the fpq. [16]) For this reason, it is difficult to see how to use 
the standard, usually very simple, "inversion relation" method [4, 23] 
to obtain the free energy. 

Here we have addressed this problem. We have shown that the 
inversion and rotation relations (37) and (39) give (41), and that this 
is sufficient to deduce the "infinite lattice functional relations" (22) - 
(28), which from this point of view are simply definitions and direct 
consequences of the definitions. 

We then made two assumptions as to the analyticity properties of 
the functions, and showed that these were sufficient to evaluate the free 
energy of the model by Wicner-Hopf factorization. The first assump- 
tion is very plausible and perfectly typical of the inversion relation 
method, which needs such assumptions to complete it. The second is 
by no means so obvious, implying as it does the symmetry (49). We 
have presented an argument for it - namely that Mq is not a branch 
cut in the tg plane for the function Tj^{tq) = T2{tg) • • • T2{io^~'^tq), and 
the easiest way to ensure this is to require that T2(tg) only have the 
cut Mjv-i. However, we must admit it is unlikely that we should have 
made such an assumption if we had not already known that it was 
implied by the exact finite-size functional relation method used in [7] 
and [8]. 

In Appendix B we present for reference the solution for the N = 2 
Ising case, where the standard elliptic function methods work immedi- 
ately. However, in this case Assumption 2 is redundant, being a direct 
consequence of Ti\f{tq) not having Mq as a branch cut. 

In later papers we intend to obtain the full Riemann surface, and 
all poles and zeros, of T{xp,yp\xq,yg) from (70), (71) and (75). This 
gives a very full description of the function and should greatly extend 
our understanding of it. 



Appendix A: checking the equations 

When performing the intricate calculations necessary to derive equa- 
tions such as (54), (62), one encounters many irritating and uninter- 
esting factors that are independent of q. Two checks on these can 
be obtained from the invariances discussed near the end of section 2: 
these can be used to immediately check that the powers of k and yp 
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axe correct. In particular, T2{tq) is proportional to ^, while ^g, S{Xq) 

are proportional to Up^^^ The functions T(xq,yq) and T[yq,Xq) 
do not contain any such external i/p factors. 

We are still left with constant factors coming from products of 
powers of io;-' and 1 — uj^ . To check these, the author has found it 
helpful to consider the low-temperature limit in which k' 0, Xp and 
\q are finite, and Xp,yp,Xq,yq — > k~^^^ , tp,tq — > k~^^^ . One needs 
to first eliminate zero differences like yq — Xp, tp — tq by using exact 
formulae such as 

— = ^^(-*- ~ W) _^N ^ k'{Xp - Ag)(l - XpXq) 

(Al) 

One can also use formulae such as 

N-l N-1 

(though we have tried to arrange the expressions so that only the first 

of these is needed) . 

In this limit T{xq,yq) = T{yq,Xg) = 1 and 

T2{J-Hq) = (1 - a;^)/(A:2/^y^^ for j = 1, . . . , iV - 1 , (A2) 
while from (23) and (9), 

/ -1 , X — XpXq 

""^^'^ '-^^ - Nk^y^pX, 



2 

p/^q 



(A3) 



It follows from (28) that 

iV-l 



fe'(l — XpXq) 

ky^xl^ 



p 

Six,) = C.(l - A,A,)-- = I ^^V } (A5) 



and from (57) that 

kyp Xp 



N-l 



The rhs of (58) is now 

^ 2N-2 

k'{Xp — Ag)(l — ApAqr) 

^UpXp Xq 

and we see that (58) is indeed satisfied in this low-temperature limit. 
Also, from (60), 

Apq = Bpq = Q , (A6) 

so (A5) agrees with (62). 
From (65) and (66), 

Dpq = l , Ppq = (Ap/A,)(^-i)/2 , (A7) 

so (67) is satisfied. 
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Appendix B: The Ising case 



When N = 2, as in [15], the parameters aq,bq,Cq,dq of ([18]) can be 
written as 

aq,bq,Cq,dq = -H (Uq) , -Hi{Uq) , @i{Uq) , Q{Uq) , (Bl) 

from which it follows that Xq,yq,iiq are 

Xq = -Vksnuq , yq = -Vk ^^^'^ , Uq = ^ . (B2) 

an.Uq an.Uq 

Here H{u),Hi{u),Q{u),Qi{u),snUq,cnUq,AnUq are the usual Ja- 
cobi elliptic functions of modulus k. For [xp, yp), (xg, (yq) G Dq we can 
restrict Up and Uy to lie in the corresponding rectangle in Figure 4, 
i.e. 3K/2 < ^{uq) < 5K/2, -K' < Q{uq) < K' . (The variables and 
functions Xq, yq, tq, Xq, T{xq, yq), T2{tq), S{Xq) in this Appendix are all 
periodic of period 2iK' .) Similarly, {xq,yq) G "Dg (or D'l) if Uq lies in 
the rectangle Pq (or T)[) in Figure 4. The vertical lines in between 
the rectangles are where lA^I = 1. 

Defining 

sn u 

scdu= , hAu) = Hi(u/2)ei(u/2) , 

cn u dn n 

it follows that 

Wpq{l) = k'scd (^^y^) , Wpq{l) = k'scd (|) , (B3) 

where K [K') is the complete elliptic integral of modulus k {k') and 

u = Uq — Up (B4) 

and by definition VFpg(O) = TFpg(O) = 1. 

We note that the Boltzmann weights depend on Up, Uq only via 
their difference Uq — Up. This is the "difference property" for this 
model (it only holds for N = 2). The partition function, transfer 
matrix and hence T{xq, yq) therefore also depend on p and q only via 
u = Uq — Up] in particular we can write T{xq,yq) as T{u). 

From (10), 



/ti(0 )9i(0)g i(i^) 
hi{u)^ 



rpq = riu) = j^j^ . (B5) 



Incrementing Ug by K takes Xq, yq to yq, —Xq. Hence the inversion 
relation (25) becomes 

T{u)T{u + K)= r{u) . (B6) 

Define 
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Figure 4: The complex Uq plane in the Ising case, showing the 
rectangles corresponding to Vq,V'q, V'^. 



Our assumption 1 implies that log T(?i) is analytic in the vertical strip 
< 3fi(it) < K in the complex tt-plane, and periodic of period 2iK' . 
Hence in this strip it has a Fourier expansion of the form 



oo 



logr(«) = Pm^"" ■ (B8) 

m=— oo 

Taking logs of (B6), using this expansion and equating Fourier coeffi- 
cients, we obtain 

which can be shown to be Onsager's result [1] for the free energy of 
the Ising model. If we define 

°° ^^ \4m+l (-\ „/4m+3 \4m+4 

9K.Z) - _ ^,4m+2^^4m+3 _ ^/4m+4^^4m+4 ' ^^^^> 

then we can write (B9) as 

T{u) = g{z)g{ci/z)/[g{l)g{ci)] , (Bll) 

which manifests the rotation symmetry (39), i.e. T{u) = T{K — u). 
The auxiliary functions T2{tq), S{Xq) simplify if we define 

f{u) = hi{u)hi{K-u)T{u)/[hi{0)hi{K)] , (B12) 

then 

_ 2 hl{Kf f{Uq - Up) f{5K -Uq- Up) 

'^ '^ " //iK)2GK)eiK) ' ^^'^^ 
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5^ ^ ^ 4 /ti {Kf f{u, - Up) f{4K -ug- Up) 

mu,)@{u,)f 

The first of tfiese formulae explicitly manifests the invariance of T2{tq) 
under Uq — > bK — Uq, corresponding to interchanging Pq with V'q 
and Xq with yq. The second manifests the invariance of S{\q) under 
Uq AK — Uq, corresponding to negating Xq,tq in Vq, while leaving 
liq, Xq, Uq Unchanged. Note that 

/ti(K)2 = i7i(o)2e(o)ei(o)/2 . 

There are further simplifications from using T{u) instead of T{u). 
The original inversion and rotation relations (37), (39) become 



Hi{0)' 



T{u)fi-u) = , f{u) = f{K - u) , (B15) 



while (Bll) becomes 

f{u) = e^-(^--y^^^'giz)giq'/z)/[g{l)giq')] , (B16) 
where g{z) is a considerably simpler function than g{z), being 

oo / 1 _ /2n-l \ 2n 

It is also true that 

T(,A - -^U^^ nu/4K' TT (i-g Z) [i-q z ) 

^""^ H.iO)'' |i^(l-g'2n-l^-l)2n-l(l-g'2n^)2n "^^^^^ 
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